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ABSTRACT 

The ISO/LWS far-infrared spectrum of the ultraluminous galaxy Mkn 231 shows OH and H2O lines 
in absorption from energy levels up to 300 K above the ground state, and emission in the [O I] 63 fiin 
and [C II] 158 /xm lines. Our analysis shows that OH and H2O are radiatively pumped by the far- 
infrared continuum emission of the galaxy. The absorptions in the high-excitation lines require high 
far- infrared radiation densities, allowing us to constrain the properties of the underlying continuum 
source. The bulk of the far- infrared continuum arises from a warm (Tdust = 70—100 K), optically thick 
(rioo/mi = 1 — 2) medium of effective diameter 200-400 pc. In our best-fit model of total luminosity 
Ltr, the observed OH and H2O high- lying lines arise from a luminous (L/Ltr ~ 0.56) region with 
radius ~ 100 pc. The high surface brightness of this component suggests that its infrared emission is 
dominated by the AGN. The derived column densities iV(OH) > 10 17 cm" 2 and N(R 2 0) > 6 x 10 16 
cm -2 may indicate XDR chemistry, although significant starburst chemistry cannot be ruled out. The 
lower-lying OH, [C II] 158 fim, and [O I] 63 lines arise from a more extended (~ 350 pc) starburst 
region. We show that the [C II] deficit in Mkn 231 is compatible with a high average abundance of C + 
because of an extreme overall luminosity to gas mass ratio. Therefore, a [C II] deficit may indicate a 
significant contribution to the luminosity by an AGN, and/or by extremely efficient star formation. 

Subject headings: galaxies: abundances — galaxies: individual (Mkn 231) — galaxies: ISM — galaxies: 
starburst — infrared: galaxies — radiative transfer 



1. INTRODUCTION 

The peculiar ultraluminous infrared galaxy (ULIRG, 
L m > 10 12 L Q ) Markarian 231 (Mkn 231, 12540+5708) 
is the most luminous infrared galaxy in the local 
universe, with a 8-10 00 /im luminosity of 3.2 x 10 12 L© 
(jSanders et al.ll2003f ). and may be a representative ex- 
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ample of the l ink between AGNs and nuclear starbursts 
(|Scovilldl2004f ) . A QSO-like nucleus is evident from many 
observations: optic all y it is cla ssified as a Type 1 Seyfert 
(iBoksenberg et all Il977t ICutri. Rieke. fc Lebofskvl 
11984 iBaan. Salzer. fc Lewinterl I1998D . it exhibits UV 
through IR polari zation and broad absorption lines 
(Smit h et alj 119951). it has compact X-ray emission 
(e.g.. iGallagher et al.l I2002D and extrem ely compact 
mid-infrared emission (Soif er et alJ 120001 and in the 
radio it is variable and possesse s a parsec scale je t 
(lUlvestad. Wrobel. fc CariUl H99l iTavlor et all fl999h . 
Nevertheless, there is also evidence of a compact star- 
burst in these results as well as in VLA observations of H 
I 21 cm absorption (ICarilh, Wrobel, fc Ulyestadl 11998). 
near-infrared observations dTacconi et al.l l2002f) , and 
millimeter CO interferometry (jBrvant fc Scovillel 1 19961 : 
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iDownes fc Solomonl [19981 hereafter DS98). Estimates 
for the starburst luminosi ty range from 1/3 to 2/3 of 
the bolometric luminosity (|Davies et al,ll 20041 DS98). 

Molecular observations have provided important clues 
about the concentration and kinematics of the gas in Mkn 
231. DS98 showed the presence of an inner nuclear disk 
of radius ~ 460 pc in CO (2-1), and a more extended 
disk with lower brightness. Most of the molecular gas 
has been found to be dense (~ 10 4 cm" 3 ) and warm 
(~ 70 K) from rec ent observations of CO and HCN sub- 
millim eter lines dPapadopoulos. Isaak. fc van der Weril 
[20071 hereafter PIW07). lLahuis et alj (|2007h have irT 
ferred embedded starburst chemistry in Mkn 231 and 
other ULIRGs based on mid-IR Spitzer observations of 
ro-vibrational bands of war m/h ot HCN and C2H2, while 
iGracia-Carpio et al.l (|2006D and lAalto et all (|2007l) have 
inferred XDR chemistry and/or radiative pumping based 
on anomalous intensity ratios of millimeter lines of HCN, 
HNC, and HCO+. 

The bulk of the luminosity in ULIRGs is emitted 
at far-infrared (FIR) wavelengths, where a number of 
molecular tracers are detected, mostly in absorption. 
Prominent lines of OH and H2O were detected using 
ISO/LWS in the FIR spectrum of Arp 220, along with 
absorption features by radicals such as NH and CH, 
revealing a chemistry that may be indicative of PDRs 
with plausible contribution by shocks and hot cores 
((Gonzalez- Alfonso et al.l [2004 hereafter Paper I). How- 
ever, t hose species are also expec ted to be enhanced in 
XDRs (|Meiierink fc Spa ans 2005), so that the dominant 
chemistry in the nuclear regions of ULIRGs remains un- 
certain. In Paper I, the ISO/LWS FIR spectrum of Arp 
220 was analyzed by means of radiative transfer calcula- 
tions, which included a non-local treatment of the molec- 
ular excitation by absorption of FIR photons. Paper I 
showed that the population of high-excitation OH and 
H2O rotational levels, in evidence from absorption in 
high-lying lines, is pumped through absorption of FIR 
continuum photons, a process that requires high FIR ra- 
diation densities. The detection of these lines thus not 
only reveals the chemical and excitation conditions in 
the absorbing regions, it also sheds light on the size and 
characteristics of the underlying continuum FIR source 
in spite of the low angular resolution currently available 
at these wavelengths. 

In this paper we extend our approach of Paper I to 
the ISO/LWS FIR spectrum of Mkn 231, and show 
that this galaxy spectrum presents striking similarities 
to that of Arp 220. Specifically, strong absorption in 
the high-excitation OH and H2O lines is also seen in 
Mkn 231. Rotationally excited OH in Mkn 231 has been 
previously detected via the 2 n!/ 2 A-doublet transitions 
(|Henkel. Guesten. fc Baanl lT987). VLBI observations of 
the mega-maser OH emission at 18 cm wavelength trace 
an inner torus or disk of size ~ 10 pc around the AGN 
(|K15ckner. Baan. fc Garrettll2003D . and MERLIN obser- 
vations were able to map essentially the whole single- 
dish me ga-maser OH emissio n with angular resolution of 
w 0.3" (jRichards et al.ll2005l ). We analyze here both the 
FIR continuum emission and the high-excitation OH and 
H2O lines, as well as the [C II] 158 (im and [O I] 63 fim 
emission lines. In we present the ISO spectroscopic 
observations of Mkn 231. In [J3]we first analyze simple 



models for the FIR continuum emission from Mkn 231, 
and then examine how well those models reproduce the 
observed FIR emission and absorption lines, ^summa- 
rizes our results. We adopt a distance to Mkn 231 of 170 
Mpc (H = 75 km s" 1 Mpc" 1 and z w 0.042). 

2. OBSERVATIONS AND RESULTS 

The full 43-197 11m spectrum of M kn 231 (first shown 
and discussed by lHarvev et al 1999), was obtained with 
the LWS spectrome ter (jClegg et al.| [l996) on board ISO 
(|Kessler et al.lfl996T ). In Fig.[TJ it is compared with that 
of Arp 220 (Paper I) re-scaled to the same distance (170 
Mpc). The grating spectral resolution is ~0.3 /im in the 
43-93 fim interval (detectors SW1-SW5), and ~0.6 /im 
in the 80-197 /im interval (detectors LW1-LW5), corre- 
sponding to Av > 10 3 kms -1 . The lines are thus unre- 
solved in velocity space. The « 80" beam size ensures 
that all the FIR continuum and line emission/ absorption 
from Mkn 231 (CO size ~ 4", DS98) lie within the 
ISO/LWS aperture. 

The data (TDT numbers 5100540, 18001306, and 
60300241) were taken from the highly-processed data 
product (HPDP) dataset (called 'Uniformly processed 
LWS01 data'), and reduced using version 10.1 of the Off 
Line Processing (OLP) Pipeline system (|Swinvard et all 
1996;. We performed subsequent data processing, 
including co-addition, scaling, and baseline remov- 
ing, using the IS O Spectral Analysis Package (ISAP; 
ISturm et al.l I1998T ) and our own routines. In order to 
obtain a smooth spectrum throughout the whole LWS 
range, the flux densities given by each detector were cor- 
rected by multiplicative scale factors. Corrections were 
less than 25% except for detectors LW2 and LW3 (100- 
145 /im), for which the corrections were 30%. We thus 
attribute an uncertainty of 30% to the overall continuum 
level, as well as for the line fluxes. 

Figure [1] shows that the FIR spectra of Mkn 231 and 
Arp 2 20 are similar in key aspects (see also lFischer et all 
1999), in particular the prominent molecular absorptions 
mostly due to OH doublets (that will be referred to here- 
after as lines) and the lack of strong fine-structure line 
emission typically seen in less luminous galaxies. A closer 
inspection of the pattern of line emission/absorption in 
both sources is shown in Fig. [5J where the continuum- 
normalized spectra are compared. Of particular inter- 
est are the clear detections in both sources of the high- 
excitation OH n 3/ 2 7/2 - 5/2 84 /jm and n 3/2 9/2 - 7/2 
65 /zm lines, with lower level energies of 120 and 290 
K, respectively (see SJ3]). The 330— > 221 and 331-^220 
H2O 66-67 /im lines, both with lower levels at 195 K, 
are also detected in Mkn 231, as well as the tentatively 
identified 2 2 o^ln line at 101 /im. It is likely that the 
increased noise level at A > 160 /xm is responsible for 
the non-detection of the high-excitation H1/2 3/2 — 1/2 
OH line in Mkn 231, which is seen in strong emission in 
Arp 220. While the high-excitation OH and H2O lines 
at 65-67 /jm are of similar strength in Mkn 231 and Arp 
220, the H2O lines at longer wavelengths are undoubtly 
weaker in Mkn 231, as seen for the 322— ^n, 22o^ln 
and 221 — >lio H2O lines at 90, 102, and 108 /im, re- 
spectively. The weakness of the latter lines in Mkn 231 
suggests that the region where the high- lying H 2 lines 
are formed is relatively weak in the far-IR continuum at 
A = 90 — 108 /im. The Mkn 231 spectrum thus suggests 
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that a warm component, with relatively weak contribu- 
tion to the far-IR continuum at A > 80 /im, is responsible 
for the observed high-excitation absorptions f H3.3p . Ta- 
ble [1] lists the line fluxes, continuum flux densities at the 
corresponding wavelengths, and equivalent widths for the 
lines detected in Mkn 231. 

In the case of Arp 220, we used high-spatial resolution 
continuum measurements available in the literature to 
infer that Arp 220 is optically thick e ven in the submil- 
limet er continuum (Paper I; see also iDownes fc EckarH 
2007). The steeper decrease of the flux density with in- 
creasing wavelength in Mkn 231, however, suggests that 
it has lower FIR continuum opacities (Fig. [T]) . This ex- 
pectation is further reinforced by the detection in Mkn 
231 of the [N II] 122 /jm line, a feature not seen in Arp 
220 (Fig. [3]). Other notable differences between both 
sources are that the [O I] 63 /im line is observed in emis- 
sion in Mkn 231 but in absorption in Arp 220, and that 
the ground-state 119, 53, and 79 /im OH lines are signifi- 
cantly weaker in Mkn 231 (Fig. [2]). In modeling Arp 220, 
we were forced to invoke an absorbing "halo" to account 
for these lines; in Mkn 231, no such halo is required (§7!]). 

In the spectrum of Mkn 231, the main 119.3 fan OH 
line appears to be slightly blue-shifted relative to the 
expected position, an effect we attribute to the prox- 
imity of the line to the edge of the LW3 detector. 
There is a nearby weaker red-shifted feature, at 120 /im, 
which coincides with the expected position of the ground 
n 3 / 2 5/2 — 3/2 18 OH line, and appears as a marginal fea- 
ture in both the "up" and "down" grating scans. How- 
ever, the limited signal-to-noise ratio ((1.0 ±0.4) x 10~ 20 
W cm -2 ), the narrow appearance of the feature (« 0.42 
/im), and the fact that it is not blue-shifted as the main 
line, make that assignment only tentative. In Arp 220, 
the main OH line is not shifted because it does not fall so 
close to the edge of the detector, as a consequence of the 
lower red-shift of the source. In Arp 220, a red-shifted 
shoulder appears at 120 /im, suggesting the possibility 
that 18 OH may be responsible for it (Paper I). We can- 
not however be certain that 18 OH is detected in any of 
these sources, but given the high 16 OH column densities 
we derive in some of our models below f §3.3p and the fact 
that values of the 16 OH/ 18 OH ratios below the canonical 
value of 500 may be expected in regions where the ISM 
is highly processed by starbursts (Paper I), our tenta- 
tive identification should be followed up with future Her- 
schel Space Observatory observations with higher spec- 
tral resolution and sensitivity. Finally, the spectrum of 
Mkn 231 shows a broad feature at the position of the 
IT/2 - n 3/2 3/2-3/2 OH line (53 /mi). We note that 
the blue-shifted side of this absorption is coincident with 
the OH n 3 / 2 11/2 — 9/2 line, with a lower level energy 
of 511 K; however, the proximity of this spectral feature 
to the edge of the SW2 detector precludes any definitive 
assignment. 

The FIR detections of both NH and NH 3 in Arp 220 
were reported in Paper I. NH 3 was also detec t ed vi a 
the 25 GHz inversion transitions bv lTakano et~ail ([2005), 
who derived a NH 3 column density six times higher than 
our value. The difference likely arises because of the 
high FIR continuum opacities in Arp 220, which cause 
the observed FIR absorptions to trace only a fraction of 
the total gas column. Since there are no such extinction 
effects at 25 GHz, the NH 3 inversion transitions are ex- 



pected to trace higher NH 3 column densities. Figure [2] 
shows that, by contrast, the NH 3 lines are not detected 
in Mkn 231, although the relatively high noise at 125 /im 
does not rule out future detection of NH 3 with Herschel 
at a level similar to that of Arp 220. 

There are two marginally-detected (2.5ct level) spec- 
tral features seen at 153.0 and 152.3 /im, in the Mkn 231 
spectrum (Fig. [3]) . Although close to the expected posi- 
tion of the main NH feature at 153.22 /im, the 153.0 /im 
feature appears significantly shifted by 0.25 /im from it, 
and better coincides with the position of the OH + 2 3 — 12 
line. Also, the 152.3 /im feature lies at 0.1 /tm from the 
expected position of the OH + 22 — li line. In Paper I, we 
also suggested that OH+ could contribute to the spec- 
trum of Arp 220 for two reasons: (i) our models were 
unable to reproduce, using NH and NH 3 , the observed 
strong absorption at 102 /im, which coincides with the 
expected position of the OH + 34 — 2 3 line; (it) there was 
an absorption feature at 76.4 /im that, if real, could be 
attributed to the OH + — 3 3 transition. Since OH + has 
never been detected in the galactic interstellar medium 
or that of any galaxy, here we only highlight the intrigu- 
ing possibility of its detection in two ULIRGs. Sensitive, 
higher-resolution Herschel observations are needed to re- 
solve this tantalizing speculation. 

The luminosity of the [C II] 2 -P 3 /2 — 2 -Pi/2 fine-structure 
line at 158 /tm is 2.5 times stronger in Mkn 231 than in 
Arp 220, but given the higher FIR luminosity of this 
source (Fig [T]) , the [C II] to FIR luminosity ratios are 
rather similar, with values of 2.5 x 1 0~ 4 and 2.1 x 10~ 4 fo r 
Mkn 231 and Arp 220, respectively (|Luhman et al.l2003h . 
These are among the lowest values found in galaxies, il- 
lustrating the so-called "[C II] deficit" found in ULIRGs. 
The [C II] line emission from Mkn 231 is analyzed in §3.41 

3. ANALYSIS 

3.1. Models for the far-infrared continuum 

Figure [4] illustrates several ways that the FIR to mil- 
limeter continuum can be fit and interpreted. We first 
modeled (model A in Fig. [4ji) the far infrared source in 
Mkn 231 as an ensemble of identical dust clouds each 
of which is heated by its own single central luminosity 
source. The representative cloud is assumed to be spher- 
ical, with radius R c , and is divided into concentric shells 
whose dust temperature s are computed from the balance 
of he ating and cooling ( Gonzalez- Alfons o fc Cernichard 
1999). We used a mixture of silicate and amorphous car- 
bon g rains with optical c onstants from iPreibisch et al.l 
(1993) and iDrainel (119831). The s tellar continuum was 
taken from iLeitherer et al.l (|1999f ) . but results depend 
only weakly on this choice because the intrinsic contin- 
uum is absorbed by the dust and re-emitted at infrared 
wavelengths. Once the equilibrium temperatures are ob- 
tained for each shell, the resulting continuum emission 
from the cloud is computed, and multiplied by N c , the 
number of clouds in the source required to match the ab- 
solute flux densities. This scaled spectrum is shown in 
Fig. H^,. The other three models (B, C, and D, shown 
in Fig. |4p-d) use grey-bodies with uniform dust temper- 
atures Tg to characterize the continuum em ission (e.g., 
iRoche fc Chandlerill993l : lArmus et al.ll2007[ K 

Assuming that the individual clouds do not overlap 
along the line of sight, our results do not depend partic- 
ularly on the radius or luminosity adopted for the model 
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individual cloud because identical results are obtained if 
R c is multiplied by a factor of a, the luminosity by a 2 , 
N c by aT 2 , and the continuum opacity is kept constant 
(see Paper I). The models are thus characterized by the 
luminosity of the whole ensemble, the radial opacity of 
the clouds at a given wavelength (which we adopt to be 
100 /im: T"ioo^m)i and the equivalent radius of the source, 
1/2 

defined as R eq = N c R c . These parameters are listed 
in Tabled 

In model A, the individual clouds are optically thin so 
that some degree of cloud overlap would yield a similar fit 
to the continuum while decreasing the value of R eq . For 
instance, if the clouds are distributed in a spherical vol- 

1/3 

ume, R eq = N c R c giving R eq = 400 pc for clouds with 
R c = 20 pc. However, the predicted opacity through the 

1/3 

modeled region, iV c Tioo^m , will be much higher than 
that of an individual cloud, and this physical situation 
is already described in models B-C where higher opaci- 
ties along the line of sight and a more compact region of 
FIR emission are assumed. In order to avoid this model 
redundancy, we choose our continuum models such that 
an individual "cloud" describes the characteristic contin- 
uum opacity (rioo/jm in Table [2]) and dust temperature 
through the whole region (disk), so that the resulting 

1/2 

extent of the FIR emission is R eq = N c R c . 

The observed continuum can be reproduced from 
model A's cloud ensemble that is optically thin in the 
FIR. Model A also predicts that the starburst domi- 
nates the continuum for A > 15 /xm, while the torus/disk 
around the AGN would then dominate the mid-infrared 
continuum, i n quali tative agreement with the models by 
iFarrah et al.l (2003) . The equivalent radius of the star- 
burst is slightly larger than the radius of the outer disk 
observed by DS98. Because Tioo/xm is low and R eq is high, 
this model predicts that the FIR radiation density is low, 
a prediction that is not consistent with our models of the 
observed OH line strengths f ^3.2p . 

As both Tioo/xm and Td are increased in models B and 
C, the radiation density increases and, therefore, the 
equivalent size required to reproduce the observed emis- 
sion becomes smaller. As a consequence, models B and C 
predict increasing compactness of the dust clouds respon- 
sible for the FIR emission, with R eq = 400 and 200 pc 
respectively. With a single-component model, however, 
R eq cannot be reduced more than in model C without de- 
grading the quality of the fit. However, a two-component 
model as shown in D is able to reproduce the FIR emis- 
sion, invoking a quite compact (~ 100 pc) and warm (100 
K) component (D waxm ), and a colder and more extended 
one that dominates at A > 80 /xm (-D C oid)- 

A convenient way to characterize the radiation den- 
sity in the modeled regions is to compute the radiation 
temperature at 100 /jm from 



T rad (100 /xm) 



hv 



k In 



1 



C 2 F 10 t m 



(1) 



where = ttN c R 2 /D 2 is the solid angle subtended by 
the modeled source, Fion^m is the predicted flux density 
at 100 /xm, and other symbols have their usual mean- 
ing. T ra d(100/im) is also listed in TableEJ together with 
the gas mass, luminosity, and fraction of the bolometric 
luminosity for each model. The calculated gas masses 



assume a gas-to-dust mass ratio of 100. In all cases, 
they are lower than the dynamical masses determined 
by DS98 when R eq is identified with the radial extent of 
the source (and therefore compatible with the inferred 
rotation velocities in the disk). Our inferred masses are 
in models B — D consistent with the mass inferred by 
PIW07, but are in all cases higher, by at least a factor 
of two, than the gas masses obtained by DS98. This 
discrepancy may be explained in at least four possible, 
different ways: (i) the physical radial extent of the cloud 
ensemble, which accounts for cloud filling, is given by 
Rt = f~ 1 ^ 2 Req, where / is the area filling factor, so that 
R eq is a lower limit of Rt', (ii) our calculated masses de- 
pend on the mass-absorption coefficient for dust, which 
we have assumed to be Ki300Aim = 0.33 cm 2 g _1 based 
on a mixture of silicate and amor phous carbon grains 
(jPreibisch et alJl!993t lDraindll98l . but could be up to 
a factor ~ 6 highe r if the dust is mainly compo sed of 
fluffy aggregates ()Kruegel fc Siebenmorg en 1994); (Hi) 
the gas-to-dust mass ratio may depart significantly from 
the standard value of 100; (iv) the masses derived by 
DS98 for Mkn 231 could be lower limits in the light of 
the submillimeter CO emission reported by PIW07. A 
combination of these factors may explain our higher val- 
ues. 

The luminosities in Table [2] account for 50-80 % of 
the observed 8 — 1000 /im infrared luminosity. Model A 
implicitly assumes that the calculated luminosity has a 
starburst origin; the luminosity from model B and from 
the cold component of model D are also attributable to 
the starburst in view of the spatial extent of the mod- 
eled source. Since model C and the warm component of 
model D are more compact, a combination of AGN and 
starburst contributions is more plausible. The surface 
brightness in model C is 4 x 10 12 L Q /kpc 2 , a factor of 
2 higher th an the peak g l obal v alue found in starburst 
galaxies by iMeurer et al.l {1997), suggesting an impor- 
tant (but uncertain) co ntribution by the AGN to the 
observed FIR emission (Soi fer et alJ [200l . Also, the 
luminosity-to-mass ratio of 500 L ^/M^ coincide s with 
the uppermost limit proposed by IScovil lc (2004) for a 
starburst. The very high surface brightness (1.3 x 10 13 
L Q /kpc 2 ) and luminosity-to-mass ratio (~ 3300 Lq/M q ) 
of the warm component of model D (D walm ), as well as 
its compactness, persuasively indicate that this compo- 
nent is most probably dominated by the AGN. The most 
plausible relative contributions by the AGN and the star- 
burst to -D W arm are discussed in SJU 

In summary, different approaches can be used to suc- 
cessfully fit the observed FIR continuum emission, with 
the properties of the clouds that emit that radiation in 
these approaches spanning a wide range of possible phys- 
ical scenarios. But ISO/LWS has provided us with spec- 
troscopic information, and we show next how the ob- 
served high excitation OH and H2O lines impose impor- 
tant constraints on these continuum models. 

3.2. Equivalent widths 

We analyze the OH equivalent widths assuming that 
the OH molecules form a screen in front of the IR source. 
The strengths of the U 3/2 7/2-5/2 and 9/2 - 7/2 OH 
doublets at 84 and 65 /im, enable us to conclude that the 
excited OH covers a substantial fraction of the FIR emis- 
sion region. Assuming that each line of the 84 /xm dou- 
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blet absorbs all the background 84 /im continuum over 
a velocity range of 250 km s _1 along each line of sight, 
and that there is no significant re-emission in the line, 
the covering factor is ~ 50%. This value may be consid- 
ered a lower limit for the following r easons. The submil- 
limeter CO line profiles shown by iPapadopoulos et alJ 
(|2007f) have FWHMs of 200-250 km s"\ and the lines 
are expected to be broadened by velocity gradients and, 
in particular, by the disk rotation; therefore, the veloc- 
ity range of 250 km s _1 assumed above is probably an 
upper limit. DS98 inferred local turbulent velocities of 
up to 60 km s _1 at inner radii (100 pc) and decreasing 
as r~ 03 . If we adopt an intrinsic Gaussian line profile 
with the highest value of the turbulent velocity, AV = 60 
km s _1 , and saturate the 84 /im line to the degree that 
an effective width 1 of 250 km s _1 is obtained for each 
component of the doublet, the derived 84 /im foreground 
opacity at line center is ~ 50, but the high column den- 
sity required for this opacity is hard to reconcile with 
that inferred from the other observed OH line strengths 
( §3.3[) . Finally, some significant re-emission in the 84 /im 
OH line is expected because the H3/2 9/2 — 7/2 OH line 
at 65 /im that originates from its upper level is detected 
in absorption. We therefore conclude that the observed 
84 /xm OH absorption is widespread, and probably cov- 
ers the bulk of the 84 /im continuum emission regions. 
On the other hand, the opacities in the high-lying 65 /im 
line should only be moderate; for reference, if we adopt 
for each component an upper limit of 150 km s _1 on the 
effective velocity interval for the absorption at each sight 
line, the minimum covering factor for this line is then 
25%. It is therefore possible that the OH responsible for 
the 65 /im absorption does not entirely coincide with that 
producing the 84 /im absorption but is only a fraction of 
the latter, consistent with its lower energy level being at 
nearly 300 K. Nevertheless, for the sake of simplicity, we 
assume in this Section that both lines arise in the same 
region -one that, on the basis of the 84 /im OH strength, 
covers the total FIR continuum. The derived OH column 
densities will be lower limits, and the inferred properties 
of the continuum source will be associated with at least 
~ 50 % of the observed FIR emission. 
The equivalent widths W are then given by 



W = 2 x 



1 - 



B„(T ex )VL 
F x 



x / (1 - exp{-T„}) dv, (2) 



where B v (T ex ) is the blackbody emission at the excita- 
tion temperature T ex of the line, f2 = nR^ q /D 2 is the 
solid angle subtended by the source, F\ is the observed 
continuum flux density at the wavelength A of the line, t v 
is the line opacity at velocity v, and the factor 2 accounts 
for the two lines that compose a doublet. The values of 
W are positive for absorption lines, and negative for lines 
observed in emission. Equation ^ applies both to op- 
tically thin and optically thick lines. For optically thin 
lines, W is proportional to the assumed column density; 
for very optically thick lines, W becomes insensitive to 
the column density and scales linearly with the turbulent 
velocity. Based on DS98, we have adopted AV = 40 km 
s -1 , which is probably accurate within 50%. 

1 The effective width is defined here as f (1— exp{— r v })dv, where 
Tv = to X exp{— (v/ AV) 2 } and tq is the line opacity at line center. 



The fractional level populations and column densities 
(and hence the opacities), and blackbody temperatures 
required to obtain W from eq. (J2J) can be estimated from 
a given column density iV(OH), and by assuming that 
the excitation temperature T ex is the same for all OH 
transitions. This assumption is certainly only an approx- 
imation. The OH levels are pumped primarily through 
absorption of FIR photons; this is quite a general model 
result when both the ground and the excited OH ro- 
tational lines (except the 163 /im one) are observed in 
absorption (see the similar conclusion in Paper I for the 
case of Arp 220, as well as the case of the H 2O lines at 
25-45 fim in Orion/IRc2. 1 Wright et alll200Ct ). When ra- 
diative excitation dominates, the level populations tend 
toward equilibrium with the radiation field within the in- 
ner cloud, and the excitation temperatures are similar for 
all transitions in those regions. However, the observed 
absorptions are produced close to the cloud boundary 
where the OH molecules are illuminated from only one 
side and T ex is less than the radiation temperature. Due 
to trapping effects, the lines with the highest opacity 
(119 and 84 /im OH lines) remain more excited than the 
other lines at inner locations of the region, but their T ex 
decrease more steeply outwards and fall below the T ex 
of other, thinner lines close to the cloud boundary. The 
detailed, non-local radiative transfer models described in 
N3.3I show that the assumption of equal T ex is only ap- 
proximately valid for the 119 and 84 /im lines, but for the 
other lines, and in particular for the 65 and 53 /im lines, 
T ex may be higher or lower depending on position in the 
region. Nevertheless, the opacities of the 119, 84, 65, 
and 53 /im lines are mostly determined by the T ex of the 
119 and 84 /im lines. As a result, the following analysis 
of the radiation density (or, equivalently, R eq , required 
to account for the observed absorptions) is accurate, at 
least to a first approximation. 

The three parameters now required to estimate W from 
equation @ are then A(OH), R eq , and T ex . Figure [5] 
shows the expected values of W for models A, B, and 
C as a function of T ex for the 119, 84, 65, and 53 /im 
lines, and compares them with the observed values. Each 
panel assumes a value of R eq that corresponds to the 
continuum models A, B, and C described in §3.11 (see 
also Table [2]). The values of iV(OH) used for each model 
are just reference values discussed below. 

The spectral line analysis of Model A (Fig. [5^) shows 
that this scenario can be ruled out as the main source of 
FIR radiation from Mkn 231 (Fig. [4^): the continuum 
model predicts low T ra d (between 20 and 34 K for the 
different lines) but the observed absorption in the 65 /im 
line, with lower level energy of 290 K, requires a much 
higher T ex . Furthermore, the 119 and 84 /im lines are ex- 
pected to be in emission (W < 0) as soon as T ex becomes 
higher than 20 and 25 K (which, on the other hand, is 
not possible for radiative excitation). With such a low 
value of T ex , W^^m is negligible even with our adopted 
screen 7V(OH) = 10 18 cm~ 2 , a value that overestimates 
W^s^im by more than a factor of two. 

The above problems still remain to some extent in 
model B, when R eq is reduced to 400 pc (Fig. \Ejp). Here 
the dust radiation temperatures allow the lines to be seen 
in absorption up to T ex m 40 K, however is un- 

derestimated by more than a factor of 2 for the adopted 
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AT(OH) = 3 x 10 17 cm -2 , yet this column density still 
overestimates the absorption of the 53 fim line. Although 
model B cannot account for the 65 /im line strength, a 
region of similar size but lower iV(OH) could contribute 
to the observed absorptions of the 119, 84 and 53 /im 
lines. 

The single-component model that best accounts for the 
four observed OH lines is model C with R eq = 200 pc 
(Fig. [5b). The corresponding continuum model (Fig. \4)p), 
with Td = 74 K, also fits rather well the overall FIR con- 
tinuum emission. Significantly, our models in ^3.31 show 
that the excitation temperatures required to reproduce 
the observed equivalent widths, 40-60 K, are those com- 
puted at the cloud surface if the OH is excited by the in- 
frared emission from a blackbody at Td — 74 K. Finally, 
the dust temperature and gas mass (Table [2]) in model 
C are cons istent with the gas tempera ture and H2 mass 
derived by iPapadopoulos et ail (|2007| ) from the submil- 
limeter CO and HCN emission. They found that this 
warm gas component hosts most of the molecular mass in 
the galaxy. The H 2 column density, A^(H 2 ) ~ 1.5 x 10 24 
cm~ 2 , indicates high optical depths, as in the galactic 
Sgr B2 molecular cloud, but Mkn 231 is much warmer. 
If the column density in Mkn 231 is concentrated in a 
face-on disk of thickness H = 23 pc, as concluded by 
DS98, the expected density is n(H 2 ) ~ 2x 10 4 cm -3 , just 
the a mount needed to account f or the CO submillimeter 
lines ([Papadopoulos et "all [2007). On the other hand, if 
this warm and dense component is identified with the 
inner disk of radius 460 pc reported by DS98, the area 
filling factor is / ~ 0.2. In spite of the general agreement 
between our model C with other observations, a closer 
inspection of this model ( §3.3j) reveals some discrepancies 
with other OH and H 2 lines that suggest that a slightly 
modified scenario can better explain the overall observed 
absorption patterns. 

3.3. Models for OH and H 2 

Radiative transfer modeling of the observed OH 
an d H 2 Q lines was done using the code d escribed 
in iGonzalez- Alfonso fc Cernichard (|1997l Il999h . which 
computes the statistical equilibrium populations of a 
given molecule in spherical symmetry. Line broaden- 
ing is assumed to be caused by microturbulence. Our 
code accounts for a non-local treatment of the radia- 
tive trapping in the molecular lines and of the excita- 
tion through absorption of photons emitted by dust, as 
well as for collisional excitation. Both line and contin- 
uum opacities for photons emitted in both lines and con- 
tinuum are taken into account. Collisional rat e s were 
take n from lOffer. van Hemert. fc van Dishoeckl (|1994l ) 
and iGreen. Maluendes. fc McLean! (|1993h for OH and 
H 2 0, respectively. As we also found for Arp 220 (Pa- 
per I), the overall excitation is dominated by absorption 
of FIR continuum photons in all models. If shock condi- 
tions (high density and temperature) were assumed, only 
the absorption in the lowest-lying lines would be signifi- 
cantly affected. Once the continuum model is fixed, our 
results only depend on the molecular column densities 
and turbulent velocity (see Paper I for a fuller descrip- 
tion). 

As mentioned above ( §3.2j) . the observed absorption 
strengths are not sensitive to the amounts of OH and 
H 2 in the inner regions of the modeled regions, but only 



to the amounts of OH and H 2 that are close to the cloud 
(or disk) boundary. For this reason, we calculate two val- 
ues for the derived molecular column densities: A^ scr (X) 
denotes the column density for a shell of species X cov- 
ering the infrared source (i.e., the screen case), whereas 
N mlx (X) is the inferred column density for models where 
X and dust are evenly mixed (the mixed case). Evidently 
7V mlx will be much higher than N scr , but from our data 
there are only a few, non-definitive ways to discriminate 
between the alternatives. The 163 /im OH and 120 //m 
18 OH lines are stronger in the mixed case, but neither 
of these features is unambiguously detected. Neverthe- 
less, we do not find any strong arguments for thinking 
that OH and H 2 are only present on the surface of the 
disk, and so the N mlx values may be considered some- 
what more reliable. The abundances we derive below are 
based on this assumption; we revisit the "mixed" case 
when we discuss models for the [C II] line. 

Since model C (Fig. Hfc, Fig. Et, Table [2]) gives the 
best single-component fit to most of the OH equivalent 
widths, we first check if it can account for the observed 
OH and H 2 absorption features. Figure [6] compares the 
observed continuum-subtracted spectrum and the mod- 
eled results (dashed spectrum, mixed case) for the wave- 
length ranges where the signal-to-noise ratio is adequate. 
Table [3] lists the physical parameters obtained for this 
model. We have assumed a turbulent velocity AV of 40 
km s" 1 (g3?l DS98). The model fits satisfactorily the 
OH 119, 84, 65, and 53 /im lines, thus demonstrating 
the approximate validity of the simple method outlined 
in §37JJ The value of N scr {OH) = 10 17 cm" 2 is also the 
same as estimated from the equivalent widths. The H 2 
column densities are determined by the strengths of the 
330 ^2 2 i and 331— >2 2 o lines at 66-67 /im. 

Some features of model C, however, are inconsistent 
with the data. The possible emission in the OH 163 
/im line is not reproduced, and the absorption in the 79 
and 99 /im OH lines appears excessive. The model also 
predicts too much absorption in the 3 22 — >2n (90 /im), 
2ao-fl u (101 /im), 2 2 i — lio (108 /im), and 4 14 ^3 03 
(113 /im) H 2 lines. All these discrepancies suggest that 
the component that accounts for the absorption of the 
65-68 /im OH and H 2 lines is weaker than postulated 
in model C at wavelengths longer than 80 /im. 

These discrepancies may be resolved by invoking two 
different components for the FIR continuum emission, 
as in model D (Fig. [U Table [2]). The warm-compact 
component, responsible for the 65-68 /im OH and H 2 
lines, will produce weak absorptions in the 80-120 /im 
range as a consequence of the relatively weak continuum 
emission at these wavelengths. The more extended com- 
ponent will contribute to the observed absorptions in the 
53, 84, and 119 /im OH lines. The compactness of the 
warm component suggests that it is relatively close to the 
AGN, and thus we have assumed AV^ = 60 kms -1 for 
-Dwarm (Table [3]); this is the turbulent velocity found by 
DS98 around the rotation curve turnover radius of 75 pc. 
For the extended component (-D co id), AV = 40 kms -1 
is assumed. Figure [6] shows that a better fit to the over- 
all spectrum is indeed found with this composite model 
(grey line), with the lines in the 80-120 /im range brought 
down to levels compatible with observations. Also, the 
model predicts the 163 /mi OH line to be in emission. 
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The column densities in model D are significantly lower 
than in G and, to avoid the above mentioned H2O ab- 
sorptions at 90-120 /im, only an upper limit for the H2O 
column density has been derived for the extended compo- 
nent. On the basis of this improved fit to the molecular 
lines, we favor model D over model C. 

3.4. Models for fC II] and [O I] 

We now check if our preferred models for OH and H2O 
are consistent with the observed [C II] 158 /im and [O 
I] 63 /im lines, and in particular if the version with OH 
and dust coexistent (the mixed case) , which predicts high 
OH column densities (Table [3]), is compatible with the 
low [C II] t o FIR flux ratio seen and previously reported 
in Mkn 231 (|Luhman et al.l l2003). In Paper I we assumed 
a C + to OH abundance ratio of 100 (based on standard 
gas-phase carbon abundance, and the abundance of OH 
found in the galactic center and Sgr B2), and found a re- 
markable agreement between the predicted and observed 
[C II] line flux in Arp 220. In the present case we have 
searched for the required C + and O column densities in 
models G and D (mixed case) that fit the observed [C II] 
and [O I] line fluxes, with the additional constraint that 
the [O I] 145 /im li ne is not detected to a limit of 4x 10~ 21 
W cm (see also lLuhman et ah 2003). In model D, the 
C + abundance was assumed to be the same in both com- 
ponents. Both the [C II] and [O I] lines are expected to 
form within Ay ~ 3 mag of the surfaces of PDRs, and 
so we allowed the gas temperature to vary between 250 
K and 700 K fe.g., IKaufman et aljfl999h : the results for 
the [C II] line are not sensitive to Tk as long as Tk is suf- 
ficiently higher than the energy of the upper level, w 90 
K. In the regions directly exposed to the FUV radiation 
field, all of the carbon and oxygen are expected to be in 
their ionized and atomic forms, respectively. Therefore 
an O to C + abundance ratio of 2.1, corresponding to O 
and C + gas-phase abunda nces of 3 x 10~ 4 and 1.4 x 10~ 4 
(jSembach fc Savagel 119961), is imp o sed. U sing the col- 
lisional rates bv lLaunav fe Rouefil (|1977al |bl) for carbon 
and oxygen, we then searched for the density n(H) and 
column densities N mlx (C + ) and N mlx (O ) required to 
match the line strengths. The results for both models 
G and D are also in Fig. [6l and the calculated column 
densities are listed in Table El 

We could not find any satisfactory fit to the three 
atomic lines with model G. Since the continuum is opti- 
cally thick at A < 150 /im, extinction effects are impor- 
tant in the [O I] 63 /im line, and densities higher than 
5 x 10 3 cm~ 3 were required to account for it. On the 
other hand, high column densities were also needed to fit 
the [C II] line, and the final result is that all combina- 
tions that fitted both the [C II] and the [O I] 63 /im lines 
yielded more flux than the upper limit for the undetected 
[O I] 145 /im line. The model in Fig. [6] shows the result 
for n(H) = 8 x 10 3 cm" 3 and T k = 400 K. Only high 
densities (~ 10 5 cm -3 ) and relatively low temperatures 
(< 150 K) can approximately match both O lines, but 
then the [C II] line is underestimated by more than a 
factor of 2. 

In contrast, model D allowed us to find a more sat- 
isfactory fit to the three lines. As we argue below, the 
compact-warm component is expected to give negligible 
emission in the C + and O lines. The extended compo- 
nent -D C oid, with its moderate dust temperature {T4 — 47 



K) and continuum opacity, can explain the [O I] 63 /im 
line with a density and column density lower than in 
model G, thereby predicting a weak 145 /im line. The 
model in Fig. [6] uses for £> co id n(H) = 5 x 10 3 cm" 3 
and Tk = 400 K. It is remarkable that when we gen- 
erated the same model but excluded the effect of dust, 
we found that the [O I] 63 /im line was stronger by a 
factor of 4.2. The dust has a small effect on the other 
two [O I] 145 /im and [C II] 158 /im lines, but can 
lower substantially the [O I] 63 /im line because the ab- 
sorption of 63 /im line-emitted photons by dust grains 
(i.e., extinction) is more important than the additional 
pumping of O atoms by 63 /im dust-emitted photons. 
Therefore, the ratio of the [C II] 158 /im line to the [O 
I] 63 /im line, observed of order unity, is reproduced 
with higher densities and temperatu r es th an in previ- 
ous studies: Fig. 4 of Kaufman et al.1 (|1999h shows that 
J([OI]63/im)/J([CII] 158/im) « 1 for n(H) = 5 x 10 3 
cm~ 3 and Go ~ 10 2 , corresponding to a temperature 
some what higher than 100 K (Fig. 1 of IKaufman et all 
1999). Including dust leads to a similar intensity ratio for 
T fe = 400 K, or G > 10 3 . On the other hand, our model 
for D cold yields JQCII] 158/im)/J([OI] 145/im) rj 7.5, 
whi ch implies for the ab ove density Go ~ 10 3 (see Fig. 
6 of lLuhman et aUfeOOSl ). These results suggest that the 
effect of extinction cannot be neglected when modeling 
the [O I] 63 /im line in sources that are optically thick 
at wavelengths shorter than 100 /im. Self-absorption by 
foreground low-excitation clouds may also lower the flux 
of the [O I] 63 /i m line and and to a lesser extent that 
of the fC HI line dFisch er et al.|[l999l: IVastel et~afl 120021: 
ILuhman et al.ll2003l Paper I). 

To compare our results for line intensity ratios involv- 
ing the [ O II 63 /im line with the corresponding plots 
shown bv lLuhman et al.1 (|2003f) . one should use the ver- 
sions of our models that ignore dust, because t he in - 
tensities given in the models by IKaufman et al.l (|1999l ) 
are uncorrected for extinction. In our models that ig- 
nore dust, [CII]/[OI] ~ 0.4 and ([CII] + [OI])/I F ir ~ 
2 x 10~ 3 . Here, the FIR flux corresponds to that of com- 
ponent D r .c,w (Table J2l ). Using these values in Fig. 4 of 
ILuhman et all (|2003l )7 one obtains «(H) = a few x 10 3 
cm~ 3 and Go = a few x 10 3 , roughly consistent with our 
input model parameters. 

The C + and 0° column densities required to fit the 
lines in model D C oid are rather high. The averaged C + 
and O abundances are 2.2 x 10 -5 and 4.6 x 10 -5 rel- 
ative to H nuclei, respectively. If the C + to OH abun- 
dance ratio is 100 (Paper I), the mixed column density 
of OH is overestimated by a factor of 3. In any case, 
the iV mix (C + ) values indicate that an important frac- 
tion of C nuclei is in ionized form, > 15% in model 
- Pcoid for a gas phase car bon abundance of 1.4 x 10~ 4 
(|Sembach fc Savage! 1 1996[ K The corresponding mass of 
atomic gas is 8 x 10 8 M , at least one order of mag- 
nitude hig her than the amount estimated in M82 and 
NGC 278 dWolfire et alJfl99?l iKauTman et al.l[j"999h . As 
mentioned above, C + is expected to be abundant only in 
regions within Ay < 2 mag from the surfaces of PDRs, 
and the relatively high fraction of C + derived above indi- 
cates that those regions must be widespread in the Mkn 
231 disk. On the other hand, one possible interpretation 
of the high H2 column densities required to generate the 
radiation field needed to account for the OH lines, is that 
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we are viewing a number of PDRs that overlap along the 
line of sight (see also [JH) . 

We now check the possibility that a significant frac- 
tion of the [C II] 158 /im emission arises from a more 
extended region, like the outer disk found by DS98 in 
CO. For an outer disk mass of 2 x 10 9 M and den- 
sity of 10 2 cm -3 (DS98), and using an upper limit for 
the C + abu ndance of < 10~ 4 and collisional rates at 
250 K Cfrom lLaunav fc Rouefl[l977aL we calculate that 
< 1/4 of the observed [C II] line arises from this extended 
region. Therefore we do not expect much contribution 
from it, unless the C + is concentrated in clumps of much 
higher density. 

Our model results can be used to re-examine the is- 
sue of the [C II] line deficit found in ULIRGs. Con- 
trary to what might be expected from the low [C II] to 
FIR luminosity ratio in Mkn 231, relatively high aver- 
age abundances of C + are required in Z? co id to reproduce 
the [C II] line. Explanations based on self-absorption 
effects, extinction by dust, effects of dust-bounded HII 
regions, high Go/n, and high-formation rate of CO at 
the expense o f C + have been invoked to account for the 
JC III deficit dLuhman et a.l.lll998l [20031: iMalh otra et all 
120011 ; iPapadopoulos et al.ll2007l Paper I, and references 
therein). In our models for Mkn 231, the high abundance 
of C + yields the observed low [C II] to FIR luminosity 
ratio because of the high luminosity-to-gas-mass ratio 
(L/M) of our continuum models. In both models G and 
D, an overall L/M w 500 L Q /M Q is derived (Table [2]), 
where the gas mass includes both the atomic and molec- 
ular components. We argue that, if this extreme L/M 
ratio applies to the whole region from which the FIR 
emission arises, as appears to be the case for Mkn 231, 
the [C II] line strength relative to the infrared emission 
can never attain values in excess of 0.1%. The [C II] 
luminosity emitted per unit of total gas mass is given by 
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>[CII] 



1 



f- exp{E u /kT ex } 



xn c +A u ihv, (3) 



where gi and g u are the degeneracies of the lower and up- 
per level, E u is the upper level energy, A u i is the Einstein 
coefficient for spontaneous emission, and Uq+ denotes the 
number of C + ions per unit of total gas mass (including 
all components). It is assumed in eq. (J3J) that the line is 
optically thin. In the limit T ex 3> E u /k one obtains 
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where x(C + ) denotes the average C + abundance relative 
to total H nuclei. Therefore, the expected [C II]-to-IR 



emission is 



<ft[CII] 

L/M 



10-3 X ^ 



io- 



(5) 



for L/M — 500 L /M . Hence, even assuming the 
highest nearly possible x(C + ) ~ 10 -4 , the [C II]-to-IR 
ratio is not expected to be above 10 -3 for such high val- 
ues of L/M . The reason is simply the limited reservoir 
of gas, and hence of C + ions, per unit of output power. 

2 For simplicity, the luminosity accounts here for the total in- 
frared emission as derived in our models, whereas in previous stud- 
ies the FIR emission was define d between 40 and 120 /jm ; the cor- 
rection factor is less than 2 (sec Sanders & Mirabel 1996). 



High L/M values are in Mkn 231 primarily due to the 
contribution of the AGN, but could also be associated 
with high sta r formation efficiencies (maybe extreme, see 
lScovillell2004ri : if we attribute the total infrared luminos- 
ity to FUV heated grains in PDRs and use the units 
of Go (i.e. the Habing Field of 1.6 x 10~ 3 ergs cm~ 2 
s _1 ) to calculate the L/M for a PDR column density of 
A v = 10 mag with the standard gas-to-extinction con- 
version 7V(H + 2H 2 )/A„ ~ 2.5 x 10 21 cm~ 2 mag -1 (e.g., 
lDickmanl ll978). we find that the value of Go that cor- 
responds to L/M = 500 L /M is 2.5 x 10 4 , which we 
express for the general case as 







^ = 500 X (2^10 



10 mag 

A v 



& (6) 

so that high L/M values would indicate high Go val- 
ues. For PDRs with G ~ 2.5 x 10 4 and A v ~ 10 mag, 
the average C + abundance is ~ 2 x 10~ 5 (i.e., carbon 
is ionized within ~ 2 mag of the surfaces of PDRs) and 
4>[cu] / \L / M) will attain the lowest values observed in 
ULIRGs. 

In our favored model D, the component D W arm 
(L/M = 3270 L /M ), which emits most of the IR lumi- 
nosity, could only account for the [C II] line if essentially 
all gas-phase carbon were in ionized form. However, 
this seems implausible, as the C + abundance in XDRs 
is x(C+) ~ a fewxlO" 5 (jMeiierink fe Spaansll2005h . and 
FUV photons from stars are absorbed in thin shells of 
dust on the surfaces of PDRs. With the assumed col- 
umn densities in Table [31 the calculated contribution 
from -D warm to the observed [C II] emission is only 8%, 
the [O I] 145 fiuY flux is 15% of the upper limit, and 
the [O I] 63 /im is predicted in weak absorption against 
the continuum. (In this calculation, Tk = 10 3 K and 
n = 2 X 10 3 cm~ 3 are assumed, so that the atomic re- 
gion is in pressure equilibrium with the cooler molecular 
gas with T k ~ 10 2 K and n = 2 x 10 4 cm" 3 .) It is 
therefore expected that the strongest source of infrared 
radiation (the warm component) gives low [C II] and 
[O I] line emission even if one assumes relatively high 
averaged fract i onal a bundances for C + and O . As in 
iLuhman et all ([20031) . we then invoke a component that 
emits strongly in the FIR but yields we ak emission in 
the [C II] and [O I] lines. In contrast to ILuhman et alj 
( 2003), however, we leave open the possibility that -D W arm 
has a relatively high C + abundance, as expected from 
both XDR and PDR chemistry. In starburst and normal 
galaxies, the moderate star formation efficiencies (e.g., 
~15 L /M in NGC 1068: iPlanesas. Scoville. fe Mversl 
119911 ). implying lower values of Go, will not place such 
a strong limit on the C + cooling rate relative to the in- 
frared emission. We suggest, then, that the [C II] 158 
/im line deficit found in ULIRGs may be associated with 
a high L/M ratio that, allowing a relatively high C + 
abundance, indicates an important contribution to the 
FIR emission by an AGN, and/or a high star-formation 
efficiency. 

4. DISCUSSION AND CONCLUSIONS 

Simultaneous modeling of the FIR continuum and line 
absorptions is the only consistent way to derive from 
FIR observations the physical and chemical properties 
of bright-infrared galaxies, where due to high radiation 
densities, the OH and H2O lines are pumped through 
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absorption of continuum photons. In Mkn 231, the high- 
lying 65-68 /im OH and H2O lines indicate the presence 
of a compact (~ 100 pc) and warm (~ 100 K) region, pre- 
sumably located around the AGN, with high continuum 
opacity (rioo^m ~ 1). The size of this region is compara- 
ble to (though somewhat larger than) the inner nuclear 
disk where the OH mega-masers have been observed from 
VLBI observations (Klockn er et al.1 12003). and is similar 
to the spatial extent of the regions where the bulk of 
the OH mega-maser emiss ion is generated (reg ions C, 
NE, and SW in Fig. 5 of IRichards et al.ll2005f) . This 
strongly suggests that we are obser ving the FIR counter - 
part of the 18 c m OH emission (e.g.. lSkinner et aOl997f l. 
IRichards et all (|2005l) estimate iV(OH) ~ 1.5 x 10 45 
cm~ 2 in 1 pc clouds where Ay w 1.7 kms ; there- 
fore, N(OR)/AV - 10 15 cm- 2 /( km s _1 ). In both the 
warm and cold components of model D, we have ob- 
tained N scr (OU)/AV - 1.5 x 10 15 cm ~ 2 /(kms x ), in 
rough agreement with the estimates by IRichards et al.l 
(2005). The "screen" case is most appropriate for com- 
parison here because the OH megamasers amplify the 
background radio continuum emission. 

The nuclear region is surrounded by a more extended 
(R eq ~ 350 pc) and colder (~ 45 K) region, which pro- 
duces the bulk of the [C II] and [O I] line emissions and 
contributes to the low-lying OH lines. Allowing for an 
area filling factor of 0.6-0.7, the total extent of the FIR 
emiss ion is then comparable t o the size of the CO (inner) 
disk (jBrvant fc Scovilldll996l DS98). The inferred den- 
sities (n(H 2 ) ~ 10 4 cm -3 ) may explain the submillime - 
ter CO emission detected bv lPapadopoulos et alj {2007), 
and somewhat lower densities (a few x 10 3 cm -3 ) can 
account for the [O I] line emission. We do not however 
rule out the possibility that some of the FIR emission 
arises from a more extended region, as the outermost 
disk observed by DS98 in the millimeter CO lines or 
the much more extend ed region observed in soft X-rays 
(|Gallagher et al.ll2002T) . 

The FIR emission from D W arm is probably dominated 
by the AGN, as pointed out in g§J] IDavies et afl (| 2004( 1 
estimated that, within a radius of 330 pc, the starburst 
luminosity accounts for (0.8 — 1.3) x 10 12 L ; since the 
luminosity for the component -D C oid is ~ 7.5 x 10 11 Lq, 
the starburst luminosity co ntribution from D war m is most 
probably < 6 x 10 11 L Q . IDavies et all (| 2004) also es- 
timated that, within 100 pc of the AGN (the size of 
Avarm), the starburst contribution is (3.2— 4.7)xl0 11 L Q , 
in agreement with the afore mentioned upper limit. On 
the other hand, the extended (360 x 260 pc 2 ) rad i o con - 
tinuum 1.4 GHz emission detected bv lCarilli et"aH (|1998l ) 
is compatible with both the whole FIR emission arising 
from the starbur st (yielding an IR- to-radio flux density 
ratio of Q ~ 2.5. ICarilli et al.lli998l ). and only the D co id 
FIR emission arising the starburst (yielding Q = 2.1). If 
one now assumes that the starburst contribution to the 
FIR emission of D waim is 3 x 10 11 L Q , the resulting Q is 
2.3 (in agreement with the media n value of Q = 2.3±0.2 
in normal galaxies; Condon 1992, and references therein). 
Also, the surface brightness and luminosity-to-mass ra- 
tio of -Dwarm due to only the starburst component would 
then be 2 x 10 12 L Q kpc -2 and 550 Lq/M q , which are 
probably upper limits for a starburst. In summary, our 
preferred model D favors that ~ 2/3 of the bolometric 
luminosity from Mkn 231 is due to the AGN. 



Although the chemistry of the extended starburst 
-Dcoid component, where the lower-lying OH lines and the 
bulk of the [C II] 158 /im and [O I] 63 /im lines are gen- 
erated, appears to be dominated by PDRs, the dominant 
chemistry of the warm component may be of a different 
nature. Since the compactness of -D war m suggests that its 
FIR emission is dominated by the AGN, and given the 
high X-ray intrinsic lumino sity of the AGN (~ 10 44 erg 
s" 1 ) in the 2-10 keV band (jBraito et al J 12004( 1 . we con- 
sider the possibility that the derived OH and H2O col- 
umn densities can be better explained in terms o f XDR 
chemistry. Calculations by iMeiierink fc Spaansl (2005) 
show that, in dense XDRs with Nh ~ 10 24 cm~ 2 , the 
OH and H2O column densities are expected to attain val- 
ues of ~ 10 18 cm -2 , similar to the values that we have 
derived in the mixed case. Furthermore, the Compton- 
thick screen with Nh ~ 2 x 10 24 cm" 2 required to block 
the primary X-ray emission (S Braito et al.ll2004[) could be 
partially identified with -D wa rm, with similar column den- 
sity. Therefore an XDR may naturally explain the obser- 
vations of high- lying OH and H2O lines in Mkn 231. Nev- 
ertheless, both OH and H2O are also tracers of starburst 
chemistry, and a number of PDRs and hot cores overlap- 
ping along the line of sight could in principle yield the 
column densities required to explain our observations at 
least in the screen case, and could then mimic an XDR. 
We conclude that further observations are still required 
to discern between XDR and starburst chemistry. 

We suggest that the [C II] 158 /im line deficit in Mkn 
231, and also probably in Arp 220 (see Paper I), is pri- 
marily due to a component that dominates the FIR emis- 
sion but emit s weakly in th e [C I I] line, as previously 
suggested by iLuhman et alJ (|2003l ). However, we also 
suggest that this component may still be rich in C + , but 
with an extreme luminosity-to-gas-mass ratio that limits 
the [C II] luminosity per unit of luminous power in the 
continuum. Our derived luminosity to H2 mass ratio of 
500 L©/M Q is high; future studies will indicate whether 
it can be applied to other ULIRGs with similar [C II] 
deficits. 

Given the high C + column densities (a few xlO 19 
cm~ 2 ) required to account for the [C II] 158 /im line 
emission, the H2 columns derived toward the extended 
component are unlikely to be associated with non- 
overlapping star-forming regions surrounded by optically 
thick envelopes. A single PDR is not expected to have 
N(C + ) in excess of a few x 10 17 cm~ 2 . A scenario with a 
crowded population of PDRs overlapping along the line 
of sight is therefore more plausible. If a typical single 
OB stellar cluster em its up to 10 39 ergs -1 , the maxi- 
mum value found by IScoville et al.l (|2001( l in M51 and 
also the luminosity of W49, then the number of single 
clusters to produce the Mkn 231 starburst luminosity of 
10 12 L Q is > 4 x 10 6 . If these clusters are concentrated in 
a disk of radius 460 pc and thickness 25 pc (DS98), the 
mean dista nce between neigh boor clusters is only < 1.5 
pc (see also iKeto et aDll992( l . Merging of PDRs may be 
compatible with some area filling factor if, for example, 
some spiral stru cture is invoked within the disk. 

In NGC 1068. lSpinoglio et all (pOOl observed the OH 
119 /im line in emission against the continuum. On the 
basis of a possible XDR chemistry, and given the corre- 
sponding mass, density and temperatu r e that character- 
ize the nuclear region, iSpinoglio et al.l (|2005( l suggested 
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that the OH 119 /im emission line could be formed in 
that region. This option also relied on the fact that rel- 
atively weak FIR emission is expected to arise from the 
circumnuclear disk of the Seyfert 2 galaxy, so that the 
OH 119 /jm line can be excited through collisions and 
emit above the continuum in these dense and warm en- 
vironments. The situation is different for Mkn 231 and 
Arp 220. In these ULIRGs, the extremely high lumi- 
nosity arising from the nuclear region, together with the 
high concentrations of gas there, make the nuclear FIR 
emission component very bright. Any possible emission 
in the 119 /xm line will be obscured by the strong absorp- 
tion, which in these objects is ~ 30 times stronger than 
the emission feature in NGC 1068. Also, the strong FIR 
radiation density pumps higher-lying OH and H 2 lev- 
els, thus producing absorptions in the higher-excitation 
lines. The OH and H 2 molecules are, therefore, poten- 
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chemical conditions of bright infrared galaxies. 
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High-excitation OH and H2O lines in Mkn 231 



TABLE 1 

Line fluxes, continuum flux densities, and equivalent widths for the lines detected in the 

ISO/LWS spectrum of Mkn 231 



Species 


Transition 


\ a 
^rcst 


Line flux b 


Continuum c 


W d 






(/im) 


(10-20 W cm" 2 ) 


(10" 19 W cm" 2 ^tm- 1 ) 


(10- 2 /im) 


OH 


Hi/2 " n 3/2 3/2 - 3/2 


53.3 


-7.0 ±2.0 


33.5 


2.1 ±0.6 


[OI] 


3 Pi - 3 P2 


63.2 


2.6 ±0.8 


25.0 


-1.0 ±0.3 


OH 


n 3 / 2 9/2 - 7/2 


65.2 


-4.1 ±0.6 


23.8 


1.7±0.3 


H a O 


330 — 221 


66.4 


-3.3 ±0.9 


23.1 


1.4 ±0.4 


H a O 


331 — 220 


67.1 


-2.2 ±0.9 


22.7 


1.0 ±0.4 


OH 


n 3/2 7/2 - 5/2 


84.5 


-8.6 ±0.7 


13.8 


6.2 ±0.5 


H a O 


220 — in 


101 


-1.3 ±0.3 


8.70 


1.5 ±0.4 


OH 


n 3/2 5/2 - 3/2 
2 2 "2 Pl 


119.3 


-3.6 ±0.3 


5.22 


6.9 ±0.6 


[Nil] 


121.8 


1.5 ±0.2 


4.72 


-3.2 ±0.4 


[CII] 


P3/2 — Pl/2 


157.7 


3.7±0.1 


1.73 


-21.4 ±0.6 



a For OH doublets, an average for the two components is given. b Errors do not include calibration uncer- 
tainties in the continuum level. Negative (positive) values indicate that the line is detected in absorption 
(emission). Uncertainties in the continuum level are less than 30%. d Equivalent widths are positive (neg- 
ative) for lines observed in absorption (emission). 



TABLE 2 
Models for the continuum emission 



Model 


Component 


TlOOum 


Req 


JV(H 2 ) a 


T d 




M 


L 


L/L m 








(PC) 


(cm - 2) 


(K) 


(K) 


(M Q ) 


(L©) 




A 




6.8 x 10^2 


1830 


4.6 x 10 22 


variable 15 


23 


1.0 x 10 1() 


2.4 x 10 12 


0.72 


B 




0.55 


400 


3.7 x 10 23 


55 


44 


4.0 x 10 9 


1.7 x 10 12 


0.50 


C 




2.2 


200 


1.5 x 10 24 


74 


71 


4.0 x 10 9 


2.0 x 10 12 


0.61 


D 


warm 


1.1 


106 


7.2 x 10 23 


100 


85 


5.5 x 10 8 


1.8 x 10 12 


0.56 


D 


cold 


0.97 


350 


6.5 x 10 23 


47 


42 


5.3 x 10 9 


7.4 x 10 11 


0.22 



a A gas-to-dust mass ratio of 100 is assumed, together with a mass-absorption coefficient for dust of 44 cm 2 
g- 1 at 100 /im. b The temperature is calculated from the balance between heating and cooling; see text for 
details. 
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TABLE 3 

Models for molecules and atoms 



Models 




C 


D 

(warm) 


D 
(cold) 


NCHo) (cm~ 3 ) 




1.5 x 10 24 


7.2 x 10 23 


6.5 x 10 23 


AV (km s" 1 ) 




40 


60 


40 


Af scr (OH) b (cm- 


2 ) 


1.0 x 10 17 


9.2 x 10 16 


6.1 x 10 16 


Af mix (OH) b (cm 


- 2 ) 


4.5 x 10 18 


1.4 x 10 18 


9.7 x 10 17 


x(OH) c 




3.0 x 10" 6 


2.0 x 10" 6 


1.5 x 10~ 6 


Af scr (H 2 0) (cm- 


2 ) 


5.4 x 10 16 


6.1 x 10 16 


< 5.0 x 10 15 


7V mix (H 2 0) (cm 


- 2 ) 


2.0 x 10 18 


9.6 x 10 17 


< 1.0 x 10 17 


X(H 2 0) C 




1.3 x 10~ 6 


1.3 x 10" 6 


< 1.5 x 10" 7 


Af mlx (C+) d (cm" 
Af mix (0) d (cm" : 


- 2 ) 


2.3 x 10 20 


3.1 x 10 19e 


2.8 x 10 19e 


! ) 


4.9 x 10 20 


6.7 x 10 19c 


6.0 x 10 19 ° 



a Models correspond to the continuum models listed in Table [2] 
and shown in Fig. QJ; and d. b N BCI (X) is the column density 
of a shell of species X surrounding the continuum source, and 
JV mlx (X) is the column density for the case that species X is 
evenly mixed with the dust. c Abundances are given relative to 
H 2 using the A fmlx (A) values. d The column densities of C + and 
O are computed by assuming gas at 400 K, and a O to C+ 
abundance ratio of 2.1.° Equal averaged abundances are assumed 
in the warm and cold components of model D. 
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Fig. 1. — Comparison between the FIR emission from Mkn 231 (black) and Arp 220 (grey). The spectrum of Arp 220 has been re-scaled 
to the distance of Mkn 231 (170 Mpc). The position of the OH lines discussed in the text are indicated and labeled with their wavelengths. 
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Fig. 2. — Comparison between the continuum-normalized spectra of Mkn 231 (upper spectra, solid lines) and Arp 220 (lower spectra, 
grey lines). The positions of lines discussed in the text are indicated and labeled. 




Fig. 3. — Comparison between the spectra of Mkn 231 (upper spectrum, histogram) and Arp 220 (lower spectrum, solid line) around 
152 (im. The spectrum of Arp 220 has been re-scaled to the distance of Mkn 231 (170 Mpc). The position of the NH, CH, and OH+ lines 
arc indicated and labeled. 
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Fig. 4. — Different models for the FIR continuum emission from Mkn 231. The whol e set of parameters that characterize the models A, 
B, C, and D, are listed in Table [2] Flux densities at 800 and 1100 are ta ken from IIRo chc & Chandler 19931 correc ted for non-thermal 
emis sion), a t 450 /im from Rigopoulou, Lawre nce, fc Row an-Robinson (1996), at 350 /im from lYang fc Phillips! <2007D . and at A < 25 /im 
from Ricki l(T976l ). 




20 25 30 30 40 50 40 60 

Excitation temperature (K) 

Fig. 5. — Predicted equivalent widths (W) versus the OH excitation temperature (T ex ) for the 119 (solid curves), 84 (dotted), 65 
(dashed-dotted), and 53 (dashed) /im OH doublets. Results are shown for three assumed sizes of the continuum emission, corresponding 
to the continuum models A, B, and C in Table [2] The observed values are indicated with horizontal lines. 7Vi7(OH) is the assumed OH 
column density in units of 10 17 cm -2 . The equivalent widths are positive if the lines are in absorption; a drop in W with increasing T ex 
thus indicates that the corresponding line is turning from absorption to emission. 
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Fig. 6. — Continuum-subtracted spectrum of Mkn 231 compared with two models for the OH, H2O, [C II], and [O I] lines. The models 
are the single-component C (black-dashed line), and the composite model D (grey-solid line) whose parameters are listed in Tables l2l and 
[3] and whose predicted continuum SEDs are shown in Fig. [4}>d. The lines that most contribute to the modeled spectra are labeled. 



